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Excitation Of MHD Modes With Gravitational Waves: 
A Testbed For Numerical Codes 
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We consider a gravitational wave oscillating in an initially homogeneous, magnetized fluid. The 
fluid is perfectly conducting and isentropic, and the magnetic field is initially uniform. We And ana¬ 
lytic solutions for the case in which the gravitational wave is linear and unaffected by the background 
fluid and field. Our solutions show how gravitational waves can excite three magnetohydrodynamic 
(MHD) modes in the fluid: Alfven waves, and both fast and slow magnetosonic waves. Our analytic 
solutions are particularly useful for testing numerical codes designed to treat general relativistic 
MHD in dynamical spacetimes, as we demonstrate in a companion paper. 
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I. INTRODUCTION 

Relativistic magnetohydynamics (MHD) plays an im¬ 
portant role in high-energy astrophysics. Quasars, active 
galactic nuclei, X-ray binaries, and gamma-ray bursts 
(GRBs) are thought to be powered by black holes or 
neutron stars. In most cases, relativistic plasmas and 
strong magnetic fields are believed to be involved in the 
extraction of energy from the central objects. 

In the past decade, numerical codes have been devel¬ 
oped to simulate MHD fluids in general relativity [1-6]. 
However, most of the codes to date have assumed that 
the background metric is stationary and fixed. This as¬ 
sumption is valid in many astrophysical scenarios and 
it greatly simplifies the numerical calculations. Such is 
the case, for example, for gas accretion onto Kerr black 
holes whenever the self-gravity of the fluid can be ig¬ 
nored. However, the spacetime is not stationary in many 
other cases. In the ‘collapsar’ model of GRBs [7], in rela¬ 
tivistic stars undergoing magnetic braking of differential 
rotation [8], and in the catastrophic collapse of hypermas- 
sive neutron stars [9], the spacetime is highly dynamical, 
and codes based on a fixed, stationary background metric 
are not able to simulate these scenarios. 

We have recently developed the first general relativis¬ 
tic (GR) MHD code [10] that evolves the spacetime met¬ 
ric, together with the fluid, by integrating the coupled 
Einstein-Maxwell-MHD equations (the “GR MHD equa¬ 
tions” for short) in 3-1-1 dimensions without approxima¬ 
tion [11] (hereafter. Paper I). We expect that similar 
codes will be developed by other groups in the near fu¬ 
ture. To verify our code, it is important that they pass 
some nontrivial test problems involving dynamical space- 
times. There already exists a suite of standard tests for 
GR MHD codes in Minkowski spacetime [4, 5, 12], as 
well as tests in stationary curved spacetimes [2, 4, 5]. 
However, very few nontrivial problems have been pro¬ 
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posed to test a GR MHD code in dynamical spacetimes. 
The evolution of the MHD interior and vacuum exterior 
magnetic field in the background spacetime of a spher¬ 
ical, homogeneous dust ball undergoing collapse (“mag¬ 
netized Oppenheimer-Synder collapse”) has been solved 
and provides one such test problem [13]. The purpose of 
this paper is to formulate another such test problem, one 
which is analytic and very straightforward to implement 
in a pure MHD environment. We consider a gravita¬ 
tional wave oscillating in an initially homogeneous fluid 
immersed in a uniform magnetic field. The gravitational 
wave will, in general, excite three wave modes in the fluid: 
Alfven waves, and fast and slow magnetosonic waves. As¬ 
suming that the gravitational wave is weak, we linearize 
the GR MHD equations and find analytic solutions. This 
problem has been studied previously by several authors in 
other contexts (see [14-17] and references therein). How¬ 
ever, the equations in most of those papers are written in 
the orthonormal tetrad frame [18], whereas most of the 
GR MHD codes today evolve components of the MHD 
variables in a coordinate basis. The purpose of this pa¬ 
per is to derive the linearized solution simply and present 
it in a form that can be used to compare directly with 
the numerical results. We also consider a different physi¬ 
cal scenario whose solution is more suitable for numerical 
code tests than the solutions found in previous studies of 
similar problems. We recently have performed numerical 
simulations for the test problem formulated in this paper 
with our new GR MHD code and found good agreement 
with our analytic solutions. The details of our numerical 
code and simulations are reported in Paper I. 


The structure of this paper is as follows. We briefly 
summarize the GR MHD equations in Section IL In Sec¬ 
tion HI, we derive the linearized equations. Since most 
numerical codes solve the fully nonlinear equations, it is 
important to know in what regime the linearized solu¬ 
tion applies. Therefore, we list the necessary conditions 
for the linearization procedure to be valid. We then solve 
the linearized equations in Section IV and summarize our 
results in Section V. 
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II. GENERAL RELATIVITIC MHD 
EQUATIONS 

Throughout, we adopt geometrized units in which c = 
G = 1. Greek (spacetime) indices run from 0 to 3, while 
Latin (spatial) indices run from 1 to 3. The signature of 
the metric is (—h ++). 

The law of baryon number conservation gives rise to 
the continuity equation V^(pow^) = 0 , where po is the 
rest-mass density and is the 4-velocity. The continuity 
equation can be written as 

-b 9i(p*v*) = 0 , ( 1 ) 

where p* = y/—g Pqu* , v'‘ = u^ju^, and g is the de¬ 
terminant of the metric Conservation of energy- 

momentum = 0 can be written as 

dti^TgT^) + d.iV^TK) = . ( 2 ) 

The time component of this equation gives the energy 
equation, and the spatial components give the momen¬ 
tum equations. 

The stress-energy tensor of a magnetized fluid is given 

by 


T''" = pohu^u'' + Pgf^’' + TZ , (3) 

where P is the gas pressure. The specific enthalpy h is 
related to the specific internal energy ehy h = l-be-b 
P/ po- We adopt a T-law equation of state P = (T— l)poe. 
The stress-energy tensor for the electromagnetic field is 
given by 

= , (4) 

where is the electromagnetic field tensor. In the ideal 
MHD regime, where the electric field vanishes in fluid’s 
frame, it is possible to write in terms of a magnetic 
vector field as follows: 


6 ^ 


2^4^ 



(5) 

( 6 ) 


where is the Levi-Civita tensor, and , is the 

(y) 

magnetic field measured by an observer comoving with 
the fluid. In terms of the EM held stress-energy tensor 
becomes 


TZ = + -b‘^g>^‘' - b'^b'' , (7) 


where 6 ^ = b^bf^. 

The evolution of the EM helds is governed by the 
source-free part of the Maxwell equations E)i,y,A + + 

Fu\,n = 0. This equation is equivalent to 


= 0 , 


where the dual tensor F*^’-' is dehned as 

= -^^hvk\F . (9) 

We also introduce the vector = n^F*^‘' = aF*>^*, 
where is the time-like normal vector orthogonal to 
the t=constant hypersurface which satishes = — 1 , 
and where a is the lapse function. The vector is the 
magnetic held measured by a normal observer moving 
with 4-velocity (see Paper I), and is purely spatial 
(n^H^ = 0 = H*). The vector 6 ^ is related to B^ by (see 
Paper I for a derivation) 


6 ^ = - 




( 10 ) 


where P^i, = g^^ + u^u^. 

In terms of B^, the time component of Eq. ( 8 ) gives 
rise to the magnetic constraint equation 

d^iyPfB^) = 0 , ( 11 ) 


where 7 is the determinant of the 3-dimensional spatial 
metric 77 = gij. The spatial component of Eq. ( 8 ) gives 
the magnetic induction equation 

dtiVlB^) + d,[^{v^B^ - FB^)] = 0 . ( 12 ) 


In summary, the GR MHD equations for the huid and 
electromagnetic helds are the continuity equation ( 1 ), 
the energy-momentum equations ( 2 ), the magnetic con¬ 
straint equation ( 11 ), and the induction equation ( 12 ). 
The MHD stress-energy tensor is 

= (poh -p b'^)u>^u‘' + (^ + y) 9^'' ~ 


The vector b^ is calculated from using Eq. (10). 


III. LINEARIZED EQUATIONS IN THE 
PRESENCE OF GRAVITATIONAL WAVES 


Consider a gravitational wave oscillating inside 
perfectly conducting huid immersed in a magnetic held. 
In the absence of the wave, the huid is homogeneous and 
static and the magnetic held is uniform. The spacetime 
metric can be decomposed as 

gti,. = Vtiu + . (14) 

where is the Minkowski metric and is the 

metric perturbation generated by the magnetized huid. 
We assume that the wave is linear 


\hZ < \r]Z . (15) 

We wish to ignore in our computation, so we re¬ 

quire the condition 


( 8 ) 


\hZZ « \hZ 


(16) 
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in the region of interest. 

The magnitude of can be estimated as follows. 

In the absence of gravitational waves, we set up Riemann 
normal coordinates (see Ref. [19], Section 11.6) at a given 
point inside the fluid. The metric near the point is then 
given by 


where Ra^0v is the Riemann tensor. Hence we have 

^matter ^ ^ Raul)p.)x°^x'^ (18) 

and ~ Rr^, where TZ is the magnitude of the 

Rieman tensor. From Einstein’s equation, TZ is of order 
|r^,y|. Hence, we have ~ |T^i,|r^. We consider 

the regime where which implies that 

r « 1/^^ . (19) 

Equivalently, the region we consider should be much less 
than the radius of curvature. 

Condition (16) can be written as 

|7;.| < I VIA" > (20) 

where L is the size of the region we consider (i.e., the 
computational domain). Note that this condition is more 
stringent than Eq. (19) since <C |?7^;y|. We will 

consider the situation where L is of the order of the 
wavelength of the gravitational wave. Hence we have 
^ where k is the wave number of the grav¬ 

itational wave. Let po and Po be the rest-mass density 
and pressure of the unperturbed fluid and Bq be the mag¬ 
nitude of the unperturbed magnetic field. Then the mag¬ 
nitude of T^i, is of order ~ po(l+ eo) + Co + 
where eo is the specific internal energy of the unperturbed 
fluid. Condition (16) can be written as 

S < fcAo , (21) 

where ho = is the magnitude of and 

^ = Po(l + eo) + Pq + A- . (22) 

47r 

Our linearized solution is only valid for a time much 
less than the dynamical collapse time of the unperturbed 
fluid; otherwise the background spacetime will evolve ap¬ 
preciably due to self-gravity: 

. (23) 

Equations (15), (21) and (23) are the conditions that 

must be satisfied for our solution to be valid. 

We orient the axes so that the (plane) gravitational 
wave is propagating in the z direction. In the transverse 


and traceless (TT) gauge, the metric perturbation 
takes the form 


V(C) z) 


'00 0 O' 

0 h+{t,z) hx{t,z) 0 
0 hx{t,z) —h+{t,z) 0 
.00 00 . 


(24) 


where /i+ and hx are the plus- and cross-polarization of 
the gravitational wave, respectively. The inverse of 
is , where = rj^°‘r]'^^hap- The metric 

determinant is (/ = — 1 -I- 0{\h^^\^). In the presence of 
the gravitational wave, we write 

=Bl + 6B\ P = Po -f (5P , = Sv^ , (25) 

where the subscript 0 denotes the unperturbed quantity, 
and where the perturbation in density Sp will be related 
to 6P below. 

In order to expand the stress-energy tensor of Eq. (13) 
to linear order, we write 


vAr 6* 

= gy,^B'‘u^ = B^bvi , 

(26) 

vAr 6* 

= ^ + ^b*P = Bl + 5B^ , 

(27) 

47r6^ 

= Bl + 2Bo^5B^ 

+/r+[(PoA - (AA] + 2/ixSo^BoA 

(28) 

V^bx 

= B^ + 5B^ + h+B^ + hxBl , 

(29) 

V^by 

= Bl + 5By - h+Bl + hxB^ , 

(30) 

^/47^6z 

= Pg + 6B^ , 

(31) 


where Bq = a/(Pq)^ -I- (Pq)^ + (-Bq)^- The evolution of 
is determined by the linearized Einstein equation (in 
the TT gauge): 

nh+ = -Ryy) 

= 4{B^6B^ - By6By) + 16Tr 
ahx = -8n{T^y + Ty^) 

§) '<=>=» 

where we have used (13) and (26)-(31) and where □ de¬ 
notes the flat-spacetime D’Alembertian operator. The 
left hand sides of the above equations are of order k'^ho 
and hence much greater than the right hand sides by 
virtue of condition (21). Thus, we have 

Uh+ = 0 , (34) 

ahx = 0 . (35) 

To the order we are considering, the gravitational wave 

is unaffected by the magnetized fluid. 

In this approximation, the constraint equation (11) 
simplifies to 


= 4{B^6By + B^SB^) + Wtt Pq -f 




Po + ^]h+ ,(32) 


d^SB^ = 0 , 


(36) 
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and the linearized induction equation takes the form 

dtSB^ = -Bi^dj6v^ + B^djSv^ . (37) 

Since the gravitational wave is weak, the induced fluid 
motion does not develop shocks and thus the flow is adia¬ 
batic. Hence the solution to the energy equation is simply 
adiabatic compression, for which P oc pj). It follows that 
5P = (rP/po)<5po and the linearized continuity equation 
may be written as 

dt5P = -TPodjSv^ . (38) 

The momentum equation (2), setting v = i, becomes 


dt 



BoiBoj 



-d,6P - 


Bojd^6B3 


+ - 


B^djSBi 


Att 


Att 

- , (39) 


where the gravitational wave source term is 


= B§{BSd,h+ + Bld,h^)5,- 


PB-^{-Bld,h++B^d,hy,)5i« 


+{[{Blf - {B^f]d,h+ - 2B^Bp,h^ }Si^ , (40) 

and where we have used Eq. (36). Note that if Bq = 
Bq = 0, then = 0. The gravitational wave has 

no effect on the fluid in this case (in the TT gauge). 
In general, from the perturbation equations (36)-(39), 
we see that to linear order all the three-vector indices 
can be raised and lowered by the flat-spacetime metric. 
Therefore, we hereafter regard all three-vectors in the 
perturbation equations as residing in the flat spacetime. 
It is easy to show that 

= j, X Bo , (41) 

where 

Jh = - B^d.h^) , (42) 

jI = ^{B^d.h+ + Byd.h,) , (43) 

fh = (44) 


IV. ANALYTIC SOLUTION 

For simplicity, we assume that and hx take the 
form of standing waves: 

h+{t,z) = h+osinkzcoskt , (45) 

hx{t,z) = hxosinkzcoskt , (46) 

where h+o and hxo are constants. It is easy to verify 
that these forms of and hx satisfy Eqs. (34) and (35). 
One may also show that a rotation by an angle tp = 
i tan“^(/ixo/^-i-o) about the z-axis will make hx{t,z) = 
0 and h+{t,z) = ho sinkzcoskt in the new coordinate 

system, where ho = \Jh'^o + ^xo- Similarly, a rotation 

by an angle (p = —^ tan“^(/i+o/hxo) will make h+{t, z) = 
0 and hx{t,z) = ho sinkzcoskt in the new coordinate 
system. This means that setting either h+o = 0 or /ixo = 
0 will not sacrifice generality for forms of and hx 
given in Eqs. (45) and (46). However, keeping both h+o 
and hxo may still be useful for the purpose of testing a 
numerical code. 

Substituting Eqs. (45) and (46) into Eq. (40), the 
gravitational-wave source term, S'®^, is given by 

= kBo{Boh+o + Bohxo) coskzcoskt , (47) 

47rS'^^ = A:Pq(—P ph+o-f Pq/ixo) cos/czcosfct , (48) 

= -k{[{BSf-{Byf]h+o + 2BSByhxo}y< 


coskzcoskt . (49) 

The time derivative of can be written as 

dtSf^^ = -k'^coskzsinktSf^ , (50) 

= (Po"h+o + Blhxo)B^o . (51) 

= (B-hxo - Byh+o)B^o . (52) 


= -[{B-f - {Byf]h+o - ‘IB-Blhxo ■ (53) 

The solution of the linearized MHD equations is uniquely 
determined by the initial conditions. We consider the 
case where all the MHD variables take their unperturbed 
values at t = 0 and hence 


Equations (37)-(40) with the constraint (36) are the lin¬ 
earized MHD equations in the presence of the gravita¬ 
tional wave hp,i,. 

Note that our equations are expressed in a coordinate 
basis. It is easy to transform the equations to other basis 
vectors. In Ref. [16], the authors derived the MHD equa¬ 
tions in the orthonormal tetrad basis using a different 
approach. We have verified that our equations (37), (38) 
and (39) are equivalent to their equations (14), (22) and 
(24), respectively [note that the authors assume Bq = 0, 
and that a minus sign is missing in their equation (16) 
for j Since all of the GR MHD codes to date evolve 
the coordinate components of the MHD variables, we will 
continue to work in the coordinate frame for the purpose 
of setting up a code-test problem. 


5B\0,x)=0 , (5P(0,a;) = 0 , (5P(0,a:) = 0. (54) 
It follows from Eqs. (39) and (54) that 


dt[Sv - va{va ■ = 


:!GW 


E 


(55) 


where we have introduced the Alfven velocity v a-, 


Bo 


VA = 


VAjtE 


(56) 


and where 

SGW ^ gGW ^ 0) = fc coskz . (57) 
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Contracting the above equation with va gives 

0 

where we have used the fact that va ■ = 0 by virtue 

of Eq. (41). Substituting Eq. (58) to Eq. (55) gives 

aGW 

dMt^o = ■ ( 59 ) 

Taking the time derivative of Eq. (39) and using 
Eqs. (37) and (38), we obtain 

3- (“»* - + §) 

4tt 4tt 

+3^d,dk5v, + dtSf^ . (60) 

47r 

We first find a particular solution that solves Eq. (60) 
by writing 


and substitute Eq. (66) into Eq. (60), with Sf'^ = 0. 

We obtain 

— [k ■ vaY\5v + [(/« • va){k ■ Sv) — uj‘^{va • 6v)]va 

+ [(/« • va){va ■ Sv) - c^(k • 5{>)]/« = 0 . (67) 

The homogeneous solutions are the nontrivial solutions 
of Eq. (67), which are well known. The solutions consist 
of three modes corresponding to an Alfven wave and both 
fast and slow magnetosonic waves (the entropy wave is 
absent for isentropic flows). The eigenvalues and eigen¬ 
vectors when K ■ VA ^ 0 are as follows (see, e.g. [12]): 

Alfven wave (with k ■ va ^ 0)' 

ij? =u?A = (<« • va)"^ , (68) 

Sv{k) oc UA{t^) = K X Va ; (69) 

Fast and slow magnetosonic waves (with k ■ va 0)' 

up' = uJ(^ , where are the roots of the dispersion equation 

+ c^('« • va)^]uiI, + iPcKk ■ va)^ = 0 , (70) 


Sv{t,z) = Svpcoskzsinkt . (61) 

Equation (60) becomes 


[The fast (slow) mode is the one with larger (smaller) 
value of 


- (fc • VA)‘^]SVp + [{k ■ VA){k ■ Svp) - k^ivA ■ SVp)]vA 

+ [{k ■ va){va-S vp) - c^^{k ■ 5vp)]k = —S , (62) 

where k = kz. The quantity is defined as 


^2 FFp -I-i?Q/ 47r _ 2 , -.2 


— -I- c^(l — u^) , (63) 


When k-va = 0, the Alfven mode and one of the mag¬ 
netosonic modes become two linearly independent static 
modes (w^ = 0). In this case, the three independent 
modes are given as follows: 

Magnetosonic wave (with k ■ va = 0)- 


where Cg = \/rPo/(po + Po + Po^o) is the sound speed. 
Contracting Eq. (62) with k, we obtain 


2 2 2 
U = C^K , 

S'v{k) oc k ; 


(72) 

(73) 


(71) 


qz 

Si- = ^GW 

" £a-ci) 


(64) 


The other components of Svp are determined by the x 
and y components of Eq. (62). The result is 


Svj, = 




£:(1_^2)[i_(„.)2] 


^GW 




{j =x,y). (65) 


Note that this particular solution satisfies the conditions 
(li)(0, a:) = 0, but it does not satisfy the conditions (59) 
required by the initial data (54). Therefore, we need to 
add a suitable homogeneous solution (i.e., a solution of 
Eq. (60) with S'f ^ = 0) to the particular solution to 
comply with conditions (59). 

To find the homogeneous solutions, we decompose Sv 
into Fourier modes. 


Two static modes (k ■ va = 0): 

= 0 , (74) 

Sv{k) oc Va , Sv{k) oc k X Va ■ (75) 

The general solution of Eq. (60) can be written as 


Sv{t,x) = Svp{t,z)+ J d^KAA{P)UA{lP)P^'^'^ 

+ '^J d^KA„^jiK)u^j{K)F^'^'^-‘^”'A) ( 76 ) 

= SVp(t,z)F J d^K'UAif^)[CAif^)cOs{K ■ X — UlAt) 

+Da{k) sin(K • X — LOAt)] 


1 = 1 ' 


+Dmj{K) Sm(K • X — LOmjt)] 


(77) 


Sv{t,x)= J SviK)e^^^-^-^*^d^K 


( 66 ) 


where Svp{t,z) = S’VpCOskzsinkt, loa = |k • va\, 
UA{tP) = K X Va, ujmj (j =lj 2) are the two solutions 
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of Eq. (70). The vectors Uraj{n) are given by Eq. (71) 
with Wm set to u}mj- The complex functions Aa{k,) and 
Amj{i'^) are to be determined by the initial conditions. 
The functions C'/i(/«), Da{k), Cmj(,i^) and Dmj{,K,) are 
linear combinations of Aa{i^) and Amjii^)- We note that 
ua{—k,) = —ua{k) and Uraj{—K,) = Umj(,K,)- Hence 
the stationary wave solutions are obtained by setting 
Ca{-i^) = ±Ca{k), Da{-k) = ± T >^(/«), Cmj{-K) = 
and m) = . 

Since both 5vp and have coskz spatial depen¬ 

dence, we need to add to Svp a homogeneous solution 
with K = ±fc in order to satisfy the initial conditions 
[5t;(0, a;) = 0 and Eq. (59)]. Specifically, if k-VA 7 ^ 0 (we 
will consider the case k ■ va = 0 in the next subsection), 
we set 



= ^[S{k - k) + S{k + k)] , 

(78) 


= 0, 

(79) 

Crnj (^) 

= (5(«: + fe)], 

(80) 

Dmj (^) 

= 0, 

(81) 


where Aa and Amj are constants to be determined. 
Equation (77) becomes 

Sv{t,z) = Svp coskzsinkt + AAUA{k) coskzsinu}At 
+AmlUml (k) coskz sinWmit 
+Ajn 2 Um 2 {k) coskzsinu)jn 2 t ■ (82) 

We see that this solution satisfies the initial conditions 
(5i)(0, a:) = 0. The values of Aa, Ami and Am 2 are 
uniquely determined by Eq. (59). Because of the nature 
of the eigenmodes, we consider two separate cases in the 
following two subsections: k ■ va = 0 (i.e., = 0) and 

k-VA 7^0 (i.e., ^ 0). 


A. Case 1: Bq = 0 (fe ■ va = 0) 


In this case, the unperturbed magnetic field is perpen¬ 
dicular to the k vector associated with the standing wave. 
It follows from Eqs. (47)-(49), (57), (54) and (59) that 

Sv^{t,x) = Sv'^{t,x) = 0, (83) 

5t(Iu^(f,a;)|j^P = — ^^coskz. (84) 

The particular solution is given by [see Eqs. (61), (64) 
and (65)] 


Svp = z N cosfczsinfct . 


(85) 


Hence, the solution satisfying Eqs. (83) and (84) is 
obtained by adding a magnetosonic wave [Eqs. (72) 
and (73)] to 5vp-. 


Sv{t, x) = z 


Am sinkcmt + 


Shw 


sinkt 


coskz . 

( 86 ) 


gravity wave 
(standing) 



ft. 




FIG. 1: Perturbation of the magnetic field SB = {SBm + 
SBp)Bo induced by a standing gravitational wave when k ■ 
Va = 0. The unperturbed magnetic held Bo is perpendicular 
to the wave vector of the gravitational wave k = kz. Both 
stationary modes SBp and SBm have the same sinkz spatial 
dependence, but they oscillate with different amplitudes and 
at different frequencies. 


The constant Am is determined by the remaining initial 
conditions (84). The final solution is 


Sv{t, x) = —z 


S, 


GW 


f(l - c^) 


(cm sinfccmf — sinkt) cos kz . 


(87) 

The perturbations 6P and SB can be calculated by inte¬ 
grating Eqs. (37) and (38). The result is 


SP{t, x) 
5B{t, x) 


— - 2 ^ 5 —(cosfccmt — coskt) sinkz , ( 88 ) 

S{1 - cip) 


{coskcmt — coskt) sinkz , (89) 

(90) 


[6Bm{t, z) + SBp{t, z)]Bo , 


where Bq = Bq/Bq and 


SBm {t, z) 
SBp{t,z) 


S^^Bq 

W^) 


coskcmtsinkz , 


Sf^BQ 


cosktsinkz . 


(91) 

(92) 


We see that with our initial data (54), the gravita¬ 
tional wave only induces a magnetosonic wave in this 
case, but not the other two static modes. Figure 1 shows 
schematically the perturbation SB = [SBm + SBp)Bo. 
Both stationary modes SBm and SBp have sinfcz spatial 
dependence, but they have different amplitudes and os¬ 
cillate at different frequencies. 

As an interesting case, we note from Eq. (53) that the 
gravitational wave will have no effect (in the adopted TT 
gauge) on the fluid or magnetic field if SQ^r = 0, i.e. if 
Bq and B^ satisfies the equation 

[{B^q)^ - (ft")"]h+o - 2B-QBlh^Q = 0 . (93) 
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B. Case 2: Bq 7 ^ 0 (fe • va 7 ^ 0) 


In this case, we expect the gravitational wave will in¬ 
duce all the three modes. As discussed earlier, the solu¬ 
tion that satisfies Sv(0,x) = 0 can be written as 

dv{t,z) = Svpcoskzsinkt + AAUACOskzsmcuAt 
+Aml'Uml coskz sinUJmlt 
+Am2Um2COSkzSmLUm2t ■ (94) 


To satisfy condition (59), we must choose the amplitudes 
Aa, Ami and Am 2 appropriately. We note from the defi¬ 
nition of Ua, Uml and Um 2 that Ua • Uml = H-A • il'm 2 = 0 
and 

dt(uA ■ + k{uA ■ i5f3p)] coskz ,(95) 

0 ~ [Anil^mli^'^A ' f^ml) 

+ Am2‘^m2{vA ' Um2) 


+k{vA • iJ'ilp)] coskz , 

(96) 

dt(k • <5'^)]^—Q = [AmlUJml(k • Uml) 

PAm2^m2{k ' Um2) 

+k{k ■ (J'Up)] coskz . 

(97) 

Define 



(98) 

& = 

(99) 

From Eqs. (57)-(59), we have 


dt{uA ■ Sv)\^^g = ^A coskz , 

(100) 

dt(vA ■ <5v)|(^o = 0 , 

(101) 

^t(k■Sv)\^^^ = coskz . 

(102) 


Hence, matching Eqs. (95)-(97) with (100)-(102) yields 


_ ^A- k{uA ■ SVp) 


(103) 


Aynl^ml('^A * f^ml) “f Am2^m2('^A ' '^m2) 



+k{vA ■ Svp) = 0 , 

(104) 


' f*ml) 4- Am2^m2(k ' Um2) 
+k(k ■ Svp) = Cfe . 

(105) 

The last two equations determine the values of A^i and 

Am 2 - Specifically, 



Ami 

_ Cl 022 — C 2 O 12 

O 11 O 22 — O 12 O 21 

(106) 

Am2 

C 2 O 11 — Cia2i 

^11022 — ^12021 

(107) 

where 



Oil 

= A * '^ml) j 

(108) 

012 

= A ' '^m2) j 

(109) 

021 

= (A^rnli^ ' '^ml) j 

(110) 

022 

= ‘^m2(^ ' '^m2) j 

(111) 

Cl 

= -k{vA • Svp) , 

(112) 

C2 

= ^k- k{k • Svp) . 

(113) 


gravity wave 
(standing) 

'VVVVXP 


A 



FIG. 2: Gravitation-wave induced Alfven mode SBAbA when 
k ■ Va a 0- The direction of perturbation bu is perpendicular 
to both Bo and z. 



The perturbations 6P and <57?* are easily obtained by 
integrating Eqs. (38) and (37). The result is 


SP{t, z) = 


kTPo 
, Am2 


A. 


ml 


^ml(^ COSU^ml^) 


iim2(l - COSWm2t) 


tVm2 

Sv^ 

-^(1 — coskt) 


sinkz 


(114) 


dB{t,z) = k 


B^Aa 


UJA 


Ua{^ — COSLOAt) 


A. 


ml 


^ml 

Am2 

tVm2 


(Uml^O - BgUml)(l - COSUJmlt) 
Am2-^0 BQUm2)(^ COSCVm2t) 


+ -(SvpBo — B§Svp)(l — coskt) 


sinkz 

(115) 


It can be shown that in the limit B§ 0, the so¬ 
lutions (94), (114) and (115) reduce to Eqs. (87), (88) 

and (89), respectively. It is also easy to show that if 
~ GW 

k ■ S =0 [which, from Eq. (49), is the same con¬ 
dition as Eq. (93)], only the Alfven wave is excited; if 
~ GW 

(k X Va) ■ S =0, only the fast and slow magnetosonic 
waves can be excited. 

To visualize the perturbation of the magnetic field, we 
write 


SB — SBa^a 4“ SBmibmi T SBm2bm2 4“ SBpbp , (116) 

where 6 a j bmi, bm 2 i and bp are unit vectors in the di¬ 
rections of UA, uf^iBo - B^iimi, Um 2 Bo “ B^Um 2 , and 
SvpBo — BqSvp, respectively. All four unit vectors are 
perpendicular to i:. The expressions for SB a, SBmi, 
SBm 2 and SBp can be deduced from Eq. (115). Figure 2 
shows the Alfven part of the perturbation (SBAbA) and 
Fig. 3 shows the other pieces of SB. All of these four in¬ 
dependent standing modes have the sinkz spatial depen¬ 
dence, but they have different amplitudes and oscillate 
with different frequencies. 



Bo 



FIG. 3: Gravitation-wave induced magnetosonic modes 
{SBmibmi and 5Bm2bm2) and the particular solution (SBpbp) 
when k ■ va 7 ^ 0. All three unit vectors bmi, hm 2 and bp are 
perpendicular to z. All the standing modes have sinkz spatial 
dependence, but they oscillate with different amplitudes and 
at different frequencies. 

Note that for Bq = 0 there is no pressure, density 
or velocity perturbation as measured in the coordinate 
(TT) gauge adopted here. Transforming to a proper or¬ 
thonormal reference frame, where coordinate separations 
measure proper distances, the pressure and density per¬ 
turbations remain zero (they are scalar invariants) but 
the fluid acquires a familiar quadrupolar velocity pat¬ 
tern orthogonal to k [20]. Such a pattern conserves the 
comoving volume of a fluid element, and hence the den¬ 
sity and pressure. 

V. SUMMARY 

We have derived and solved the linearized equations 
for standing gravitational waves oscillating in an initially 
homogeneous, magnetized fluid (see Table I). We express 
the equations and analytic solutions in a coordinate basis 
so that the results can be compared directly with numer¬ 


ical data from a GR MHD code. 

Our linearized equations are valid when (1) the ampli¬ 
tude of the gravitational wave is small {ho <C 1), and (2) 
the magnitude of the stress-energy tensor of the unper¬ 
turbed fluid is sufficiently small that the perturbation it 
induces in the background, nearly Minkowski spacetime 
remains small compared to ho- Specifically, if the do¬ 
main in the fluid is ~ 1/fc, where k is the wave number 
of the gravitational wave, the second condition can be 
written as £ <C [Eq. (21)]. Our solution remains 
valid for a time t <C IjVS, after which the homogeneous 
background undergoes collapse. 

We consider the metric perturbations of the form 
in Eqs. (24), (45) and (46). The gravitational wave 
is a standing wave whose amplitude varies in the z- 
direction. It can be regarded as the superposition of 
two equal-amplitude, infinite traveling waves moving in 
the ±z-directions. At t = 0, the fluid is assumed to be 
homogeneous, uniformly magnetized, and at rest (i.e., 
5B = 5P = 5v = 0). The gravitational wave excites 
only a magnetosonic wave when B^ = 0. The solution 
of the perturbations is given by analytic expressions in 
Eqs. (87)-(89). When Bq ^ 0, the gravitational wave 
excites fast- and slow-magnetosonic waves and an Alfven 
wave. The solution in this case is given by Eqs. (94), 
(114) and (115). 

We have performed numerical simulations for MHD 
waves induced by linear gravitational waves with a new 
3-1-1 GR MHD code and have found good agreement with 
the analytic solutions presented here. The details of our 
numerical code, together with a comparison of numeri¬ 
cal and analytic solution for these waves, as well as the 
results of other code tests, are reported in Paper I. 
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